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Abstract 

We study the light scattering by plasmonic and dielectric symmetric trimers to investigate the 
existence of polarization-independent Fano resonances. Plasmonic hybridization theory is revealed 
to hide simple physics, and we instead provide a simplified model for hybridization to derive a 
plasmonic trimer’s eigenmodes analytically. This approach is demonstrated to accurately recreate 
full wave simulations of plasmonic trimers and their Fano resonances. We are subsequently able 
to deduce the grounds for modal interference in plasmonic trimers and the related formation of 
Fano resonances. However, by taking advantage of the generality of our simplified hybridization 
approach, we are also able to investigate the eigenmodes of all-dielectric trimers. We show that 
bianisotropic coupling channels between high-index dielectric nanoparticles are able to increase the 
capacity for Fano resonances, even at normal incidence. We finally provide the first experimental 
measurements of sharp, polarization-independent Fano resonances from a symmetric all-dielectric 
trirner, with very good agreement to the predicted response from our simplified hybridization 
theory. 
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I. INTRODUCTION 


The Fano resonance |1] has become a well-recognized interference feature in the optical 
scattering response of many nanoscale structures. It attracted considerable attention due to 
its distinctively-sharp extinction lineshape, but also due to anomalous near-field behavior for 
surface-enhanced Rahman scattering [2|, nonlinear response [3|, and enhancement of circular 
dichroism [1]. For plasmonic structures, the existence of optical Fano resonances has gen¬ 
erally been explained through the use of a hybridization theory argument [5j; an argument 
which involves subdividing a structure into two or more subsystems with known properties, 
and then deducing how their optical responses combine together. The way such optical 
responses combine, or hybridize, is regularly depicted as per molecular orbital theory: the 
modes of each subsystem are added constructively or destructively to form a bonding and 
an antibonding mode, the latter of which exhibits suppressed scattering associated with the 
Fano resonance. For the case of plasmonics, this involves treating individual nanoparticles 
as electron density distributions; a model which inherently accounts for quite comprehensive 
physics and subsequently requires nontrivial derivations for even the most simple systems, 
such as concentric spheres [6] and two-particle dimers [Tj. Using this model, even the com¬ 
paratively simple case for plasmonic hybridization of a symmetric three-particle trirner has 
yet to be found. A consequence of this complexity is that hybridization is rarely performed 
explicitly. It has instead become a conceptual explanation for experimental and numerical 
observations in even quite complicated scattering systems for which hybridization solutions 
have never been found. Here we instead begin from the premise that the dominant optical 
properties of nanoparticle geometries do not require the level of complexity found in plas¬ 
monic hybridization theory. We propose a simplified approach to hybridization when mod¬ 
eling a nanoparticle system in the dipole approximation, thereby considering only the dipole 
responses of individual nanoparticles and the coupling between them. To demonstrate the 
capacity of this simplification, we approach the outstanding problem of modal hybridization 
in symmetric trimers. In particular, the optical feature we focus on in symmetric trimers is 
the polarization-independent Fano resonance. This type of feature has previously only been 
observed in more complicated nanoparticle oligomers [8j where a conceptual hybridization 
can be performed in equivalent ways for arbitrary polarization angles. While the polarization 
independence of symmetric trimers is known M, the existence of Fano resonances mm 
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and the indeed their combination with polarization indepedence is not obvious in terms of 
hybridization, particularly because it is not possible to subdivide a trirner in a way that 
conserves its symmetry. Here we are able to resolve this problem by performing the simpli¬ 
fied hybridization procedure to find the eigenmodes of a plasmonic trirner, analytically. We 
show that the doubly-degenerate eigenmodes of trimers, those responsible for polarization 
indepedent scattering, become inherently nonorthogonal in the presence of retarded coupling 
between particles. Such eigenmodes can therefore allow polarization-independent Fano res¬ 
onances at normal incidence |14j . and we provide full wave simulations to demonstrate these 
Fano resonances in plasmonic nanoparticle trimers, with good agreement to our analytical 
eigenmode results. However, by recognizing that our simplified approach to hybridization is 
portable to other models, we further investigate all-dielectric trimers in the coupled electric 
and magnetic dipole approximation na to account for the existence of magnetic dipolar re¬ 
sponses of individual high-index dielectric nanoparticlcs. This analysis shows that previously 
reported bianisotropic coupling effects in dielectric oligomers [IB] also play an important role 
in the optical response of all-dielectric geometries at normal incidence. The result is then 
a three-fold increase in the number of doubly-degenerate eigenmodes that can be excited 
by a normally-incident plane wave in an all-dielectric trirner, when compared to the plas¬ 
monic trirner. Moreover, all these eigenmodes are nonorthogonal and subsequently lead 
to a substantial increase in the presence and magnitude of polarization-independent Fano 
resonances at normal incidence. We then verify these predictions in a radio frequency exper¬ 
iment that mimics the optical transmission through a silicon nanosphere trirner. We observe 
a number of sharp, polarization-independent Fano resonances from the all-dielectric trirner 
at normal incidence. The eigenmodes calculated from our simplified hybridization approach 
are then able to provide a quantitative recreation of these experimental measurements and 
their explicit eigenmode decompositions. In doing so, we can unambiguously demonstrate 
the validity of our simplified approach to hybridization. 


II. THE EIGENMODES OF PLASMONIC TRIMERS 

For the purpose of performing hybridization, a trirner can only be subdivided into a dimer 
and a single particle; the other option would be the trivial case of dividing it into three 
single particles. Here, we restrict ourselves to considering only the electric dipole response 
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of each individual particle, since it dominates for subwavelength plasmonic particles. We 
can therefore utilize the dipole approximation HU. where each particle’s dipole moment (p*) 
is related to the externally-applied electric field (E 0 ) as: 


P* = a^eoEotw) + a E k 2 ^ Go(r*, W 


P; 


( 1 ) 


Here, a E is the electric dipole polarizability of a particle and Go is the free space dyadic 
Green’s function, which acts on dipole moments as: 
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where k is the wavenumber, n is the unit vector pointing from r to r' and R — |r — r'|. In 
this model, the single particle and dimer eigenmodes are both known and easily deduced, 
and are shown in Fig. [l[] and [pJ- To derive how the dimer and single particle hybridize, we 
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FIG. 1. The eigenmodes of (a) a single particle and (b) a dimer with D 2 h symmetry, when assuming 
only in-plane electric dipole responses from each particle. The dimer eigenmodes are labeled 
according to their associated irreducible representation. We also show the location of the extra 
particle/s needed to form a trimer. In (c) and (d), we depict a generalized energy level diagram 
for these eigenmodes and identify the coupling channels, labeled as A-D. 

have to calculate the coupling channels between each dimer and single particle eigenmode, 
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labeled as A-D in Fig. [IJ and[l]i. To find the coupling channels, we use Eq. [I] to calculate 
the dipole moments induced in the dimer by a single particle eigenmode and project those 
onto the dimer eigenmodes (which form a complete basis for the dimer’s response), and vice 
versa. This gives us expressions for the coupling channels shown in Fig. [l/ 
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where we have normalized the single particle and dimer modes to make all channels om¬ 
nidirectional. All other coupling channels between the eigenmodes of the dimer and those 
of the single particle are zero, which is due to a geometry-induced symmetry mismatch 
between basis vectors. It is important to acknowledge that, unlike molecular orbital hy¬ 
bridization, the hybridization of scattering structures does not have a perennial set of basis 
vectors (he. atomic orbitals) with known symmetry mismatch conditions to identify recur¬ 
ring coupling/hybridization channels in different geometries. In nanostructures, the coupling 
channels have to be derived for each given geometry in order to perform hybridization cor¬ 
rectly. Importantly, our formulation of hybridization also relies on the basis vectors being 
eigenmodes of their associated subsystem. Indeed, if this was not the case, we would have to 
account for new coupling channels arising from interactions within each subsystem, which 
would substantially increase the number of coupling channels we have to consider and also 
their complexity. In any case, to now find the trimer’s hybridized eigenmodes we recognize 
that an eigenmode’s dipole moment profile must satisfy the dipole equation (Eq. [Tj) as: 

Vi = a E e 0 \vi + ^ a E k 2 G 0 (ri, rq) • v, (3) 

We can then rewrite Eq. [3] as a system of linear equations by considering the dimer and 
single particle eigenmodes [i.e. in Fig. [Tj) as orthogonal basis vectors. Moreover, we can 
consider an eigenmode, /), as a linear combination of these basis vectors: 

| Vi) = ai |p x ) + bi |B3 U ) + a |Bi g ) + di |p y ) + e* |B2 U ) + fi |A g ) (4) 
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As detailed in Appendix |Aj this basis allows us to solve Eq. [3] as a matrix equation to find 
analytical expressions for the hybridized eigenmodes of the trimer. Specifically, we obtain 
the following expressions for the hybridized eigenmodes: 
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Here all undefined coefficients (c/. Eq. [4]) are zero, and we have defined a central eigenvalue, 
Ao, for the |rq) and |tq) eigenmodes, along with a splitting function, 5, that produces their 
nondegeneracy. 
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We can categorize the hybridized eigenmodes according to their irreducible representation 
in the trimer’s D 3 h point group (see Table [T| . The |tq) eigenmode has azimuthally-oriented 
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TABLE I. Character table for the D 3 h symmetry group. The rows correspond to different irre¬ 
ducible representations and the columns are the symmetry operations. Each number is the trace 
of the associated operation’s matrix representation. [TS] 

dipole moments and the \v±) eigenmode has radially-oriented dipole moments, making 
them the trimer’s A ' 2 and A', eigenmodes, respectively p]. On the other hand, the \v%) 
and ^ 3 ) eigenmodes are the, doubly-degenerate, E' eigenmodes, which are responsible for 
polarization-independent behavior. Because we have obtained two distinct E' eigenmodes, 
which is known to be the maximum number that a ring-type oligomer can exhibit in the 
dipole approximation 0 , we know that the hybridization approach here has determined 
all of the E' doubly-degenerate eigenmodes for a plasmonic trimer. 

It is worth pausing to acknowledge an interesting result of this analysis. If we were 
to consider the doubly-degenerate eigenmodes in complex frequency space (s-space) to find 
and/or analyze their resonances [19], there are points in s-space where the expressions for \v%) 
and U 3 ) are the same, specifically: when S = 0. At such points, the two eigenmodes coalesce 
and the eigenspace of these two eigenmodes subsequently reduces in dimension, indicating a 
nontrivial topology of the eigenspace. Such coalescence points are also known as exceptional 
points and the coalescing eigenmodes are known to interchange expressions with each other, 
depending on the path taken through s-space in the vicinity of these points j2DHZ2]- In a 
trimer, the locations of four such exceptional points can be found by solving 5 = 0 for a 
complex k in Eq. [ 7 ] An investigation into eigenspace topology is well beyond the scope of 
this work and, instead, it suffices that the eigenmodes and eigenvalues in Eq. [5] are exact 
expressions for that of a symmetric trimer when considering only the electric dipole response 
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of each individual nanoparticle. Furthermore, using the expressions in Eq. [5j it follows that, 
in general, |t; 2 ) and |i? 3 ) will be nonorthogonal and written in terms of the eigenmode’s basis 
vector coefficients: 


— «2 a 3 + ^2^3 + c 2 c 3 — ^2^3 + e 2 e 3 + 72/3 (§) 

This means that polarization-independent Fano resonances could indeed be expected [Hj to 
occur in plasmonic trimers, which is interesting given neither the dimer nor single particle 
are able to exhibit Fano resonances in isolation. In other words, the presence of coupling 
channels between dimer and single particle have directly led to nonorthogonality in the 
trimer’s hybridized eigenmodes. This can be understood conceptually as the single particle 
providing a coupling channel into the |A g ) and |B lg ) basis vectors, which are orthogonal to 
the incident held and thereby act as a damped oscillators that are coupled with the remainder 
of the system. However, as is discussed in Appendix [Aj it is also important that the coupling 
channels provide retardation in propagation for Eq. [8]to not sum to zero. In any case, to now 
demonstrate Fano resonances of the nature we are describing, we can consider a symmetric 
trimer made of silver nanoparticles with a diameter of 100 nm, and we then vary the gap 
between particles from 40 nm down to touching. In Figure [2j we show the evolution of the 
extinction cross section from this trimer for a normally-incident plane wave. We provide 
two calculations for extinction: (dashed line) directly from CST Microwave Studio and 
(unbroken line) from the analytical expressions for the eigenmodes and eigenvalues in Eq. [5j 
Notably, with the exception of the quadrupole resonance at high frequencies, there is a good 
match between the eigenmode analysis and the full wave simulations. Importantly, both 
approaches observe a small Fano resonance forming as the gap between particles is reduced. 
However, the advantage of using eigenmodes is we can attribute portions of the extinction to 
individual eigenmodes. Indeed, in Fig. [ 2 ) 3 , we can see that the Fano resonance is produced 
by interference H between the two eigenmodes. Admittedly, while our model does not 
observe the single particle quadrupole response, an extension our hybridization approach to 
models that account for quadrupole responses, such as in Ref. [23], is not unfeasible. I 11 any 
case, the results in Fig. [2] are able to validate the predictions and analysis of our approach 
to hybridization for a plasmonic trimer. 
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FIG. 2. (a) Extinction cross section from a trimer made of lOOnm silver nanospheres, simulated 
using (dashed line) CST Microwave Studio and (unbroken line) the analytical eigenmodes of Eq. [5j 
Both curves show the formation of a small Fano resonance as the gap between particles, g , is 
reduced, (b) Decomposition of the extinction for the g=20nm trimer in terms of the two excited 
eigenmodes. The eigenmode profiles depict the real components of dipole moments at 400nm. 


III. THE EIGENMODES OF ALL-DIELECTRIC TRIMERS 

The analysis of the previous section corresponded to plasmonic nanoparticle trimers be¬ 
cause we assumed only electric dipolar responses from individual nanoparticles. However, 
strong Fano resonances were predicted to exist in silicon, all-dielectric, trimers [US HH- 
The hybridization of electric dipoles is not sufficient for considering all-dielectric oligomers 
because high-index dielectric nanoparticles have both electric and magnetic dipolar re¬ 
sponses [2TI - 126] . To investigate the eigenmodes of such systems, we can again employ our 
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simplified approach to hybridization, but instead use the coupled electric and magnetic 
dipole approximation [15], which will account for both electric and magnetic responses of 
individual nanoparticles. From geometric symmetry principles, we know that this approach 
will result in two distinct sets of eigenmodes for the trimer: one which will transform un¬ 
der the trimer’s symmetry operations according to the E' irreducible representation and 
the other which will transform according the E" irreducible representation (see Table |T]) . 
The intuitive distinction between the two sets of eigenmodes is that they will be excited by 
either the electric or magnetic field in a normally-incident plane wave, for E' and E", re¬ 
spectively. The different irreducible representations also means that the trimer’s symmetry 
prevents E' eigenmodes coupling to the E" eigenmodes and vice versa. If we could neglect 
the bianisotropic coupling between electric to magnetic dipoles, the E' eigenmodes would be 
those that we derived for plasmonic trimers (he. \i) 2 ) and ^ 3 ) in Eq. [5J and the E" eigen¬ 
modes would be analogous to the same eigenmodes, but constructed from magnetic dipoles 
rather than electric dipoles. However, by using the out-of-plane (z) direction, the z-oriented 
magnetic dipoles can be arranged into basis vectors that transform according to the E' irre¬ 
ducible representation and the z-oriented electric dipoles can be arranged into basis vectors 
that transform according to the E" irreducible representation, as shown in Fig. [ 3 ) 0 . Hence, 
bianisotropic coupling between electric and magnetic dipoles will allow the electric dipole E' 
basis vectors to couple with the magnetic dipole E' basis vectors, and vice versa for the E" 
basis vectors. From this point on, we shall consider only the E' eigenmodes because they are 
the eigenmodes excited by electric field and can therefore be more naturally related to the 
eigenmodes we derived for plasmonic trimers. The procedure for finding the E" eigenmodes 
is almost identical upon interchanging the electric and magnetic dipole polarizabilities. We 
begin by extending the hybridization diagram of Fig. [l]i to take into account both the new 
space of z-oriented E' responses, and the bianisotropic coupling channels. I 11 Fig. [3}i-c, we 
show the complete set of basis vectors that hybridize to form the E' and E" eigenmodes 
of an all-dielectric trimer. These basis vectors can be separated according to their even or 
odd response under the dimer’s reflection symmetry operation, because there are no cou¬ 
pling channels between the two resulting sets of basis vectors. One of these two sets is then 
sufficient to find instances of each doubly-degenerate eigenmode (i.e. being the instances 
that are even/odd under the dimer’s reflection symmetry operation). The even basis vectors 
offer a convenient opportunity to reduce the number of coupling channels because the trimer 
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Hybridization for a dielectric trimer 
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FIG. 3. The eigenmodes of (a) a single particle and (c) a dimer with D 2 h symmetry, when assuming 
both electric and magnetic dipole responses from each particle, and neglecting the z-direction. 
The (b) z-oriented basis vectors of a trimer, which become doubly-degenerate eigenmodes of the 
trimer when electric-magnetic coupling is neglected. The dimer and trimer eigenmodes are labelled 
according to their associated irreducible representation. In (d) we depict a generalized energy level 
diagram analogous to that in Fig. []}l, now with coupling channels C-F (see Eqs.[|and[§ 


basis vector is an (orthogonal) eigenmode for the dimer. The new, E and F, coupling co¬ 
efficients, shown in Fig. [3ji, can then be calculated by evaluating the magnetic and electric 
field radiated by electric and magnetic dipoles. See Appendix [B] for details. 
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We can then define an eigenmode of the dielectric trimer, |w), as a linear combination of 
single particle and dimer eigenmodes. 


w i) — a[ |p y ) + b[ |B 2u ) + c- |A g ) + d[ |E' Z ) 


( 10 ) 


Notably, the expression for the A\ hybridized eigenmode of the trimer, which has radially- 


oriented electric dipole moments, is the same as \v 4 ) in Eq. 5d This is because the eigenmode 
is unable to couple bianisotropically into any magnetic dipoles and is, therefore, unchanged 
from the previous analysis provided for plasmonic trimers. The remaining three-dimensional 
eigenspace can therefore be spanned by the two E' eigenmodes in Eq. 0 (\v 2y ) and |u 3y » 
and the |E Z ) basis vector. As such, all the remaining eigenmodes can be written as a 
linear combination of these three basis vectors, and, therefore, all remaining eigenmodes 
must transform according to the E' irreducible representation. Additionally, since the E' 
eigenmodes transform according to a different irreducible representation to that of ^ 4 ); we 
know they must be orthogonal to (tq), and we can thereby write their general form as: 
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If the bianisotropic coupling channels are negligible (E, F —> 0), the three basis vectors, 
\ v 2 y)i \ v 3 y) an d |E Z ), are the E' eigenmodes of the trimer. However, outside of this limit, we 
can use the dipole model to set up a rank 3 matrix equation to find the eigenmode solutions 
for a', d and d' (refer to Appendix B). These eigenmodes can be calculated analytically, 
but the result is cumbersome and does not provide any additional intuitive understand¬ 
ing. Therefore, it is sufficient to use a numerical approach, and this is what we do in the 
coming analysis. It is, however, worth acknowledging a conclusion here: because we have 
three linearly independent E' basis vectors, a dielectric trimer will have three E' eigenmodes, 
which one more than the plasmonic trimer. Moreover, using the analogous argument, a 
dielectric trimer additionally has another three E" eigenmodes that can be excited by the 
incident magnetic field. Therefore, the total number of, doubly-degenerate, eigenmodes 
that can be excited in a dielectric trimer with a normally-incident plane wave is six, which 
are organized into two sets of three interfering eigenmodes. For comparison, a plasmonic 
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heptamer (a central nanoparticle surrounded by a ring of six nanoparticles) has a single 
set of three interfering eigenmodes|14j. The dielectric trimer geometry therefore has double 
the propensity for Fano resonances as the plasmonic heptamer geometries. In this regard, 
the theoretical predictions for the existence of Fano resonances in all-dielectric trimers were 
originally made for silicon spherical nanoparticles mm, which exhibited Fano resonances 
in the optical frequency range. However, it is possible to create a macroscopic analogue of 
silicon nanospheres in the microwave range, specifically: using Mg 0 -Ti 02 ceramic spheres 
characterized by dielectric constant of 16 and dielectric loss factor of (1.12 — 1.17) x 1CT 4 
(measured at 9 — 12 GHz). These ceramic spheres in the microwave range therefore have very 
similar properties to silicon spheres in the optical range. As such, we are able to create a 
macroscopic analogue of a silicon nanosphere trimer. This allows us to investigate the scat¬ 
tering properties of a single trimer with much less noise than would be possible for a single 
silicon nanosphere trimer. In Fig. [IJr, we show an experimental setup. The trimer consists of 


Experimental setup 


Measurements 





FIG. 4. (a) Experimental setup for the trimer made from Mg0-Ti02 spheres and (b) the mea¬ 

surements of extinction for both the trimer and a single sphere. The trimer exhibits a pronounced, 
polarization-independent, Fano resonance at 4.8 GHz in very good agreement with (c) the simula¬ 
tion results of extinction, which were calculated using CST Microwave Studio. 

three MgO-TiC >2 spheres with 15 mm diameter and 20 mm distance between the centers of 
the spheres. To fasten together the MgO-TiC >2 particles for the experiment, we used a cus- 
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tom holder made of a styrofoam material with dielectric permittivity of 1 (in the microwave 
frequency range). To approximate plane wave excitation and receive the signal scattered to 
the forward direction, we employ a pair of, identical, rectangular linearly-polarized wideband 
horn antennas (TRIM, 1 — 18 GHz), connected to the coaxial ports of a vector network ana¬ 
lyzer (Agilent E8362C). The trimer is located in the far-held of both antennas; the distance 
from the trimer to both the receiving and the transmitting antennas is approximately 1.5 m. 
The total extinction can then be extracted from the measured complex magnitude of the 
forward scattered signal by means of the optical theorem [27]. Associated simulations of 
the experiment were also performed by using the Time Domain solver of CST Microwave 
Studio when assuming plane wave excitation on the trimer in free space. The experimentally 
measured, and numerically simulated, extinction spectra are shown in Fig. [4}r and [4] . The 
extinction spectrum of an isolated MgO-TiOo sphere was also measured and simulated. In 
Fig. [4jo, we can see a pronounced Fano resonance at 4.8 GHz that is associated with the 
suppression of extinction. By varying the orientation of the trimer, it exhibits a polarization 
independent response. The Fano resonance must subsequently be from interference between 
doubly-degenerate eigenmodes, and we can, therefore, consider the hybridized eigenmodes 
derived according to our previous analysis. For this hybridization, we define electric and 
magnetic dipole polarizabilities from the a\ and b\ scattering coefficients of Mie theory [28j. 
The Mg0-Ti02 permittivity was assumed be dispersionless with a dielectric constant of 16 
and dielectric loss factor of (1.12 — 1.17) x ICR 4 . By calculating the hybridized eigenmodes 
we get the extinction spectra shown in Fig. [5^,, which accurately reproduces the experiment’s 
extinction spectrum and Fano resonance. The discrepancy between the hybridzation theory 
and experiment at high frequencies is due to the prescence of a known [251 ESI magnetic 
quadrupole response in the individual spheres, which we simply do not take into account 
in our hybridization. Regarding, however, the Fano resonance at 4.8 GHz, we are able to 
decompose the simulated extinction spectrum into components coming from each isolated 
eigenmode as seen in Fig. [5] The extinction depicted for each eigenmode is ‘isolated’ in 
the sense that we are neglecting the extinction that can be attributed to the interference 
between eigenmodes[T4J. The decomposition clearly shows that there are multiple eigen¬ 
modes that interfere destructively to form the main Fano resonance at 4.8 GHz. Moreover, 
modal interference is occurring at very similar frequencies in both the E' and E" response; 
the Fano feature is coming from interference between EA) and |E^ 3 D, and between |E"^) 
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FIG. 5. (a) Extinction of the trimer in Fig. [4j when calculated using the hybridization procedure 
presented in this paper. We also show the eigenmode decomposition of the extinction coming from 
the E 7 and E" response of the trimer in (b) and (c). Each eigenmode’s isolated contribution to 
the extinction is shown alongside their associated dipole moment profiles, which depict the real 
components of each dipole moment at 4.8 GHz. Red arrows are electric dipole moments and blue 
arrows are magnetic dipole moments. 

and |E" 2) ). It is worthwhile emphasizing that this is a simultaneous overlap of two optical 
Fano resonances, which are symmetrically-exclusive by nature of their distinct irreducible 
representations. The likely reason for this situation is that 4.8 GHz is also the frequency of 
an individual sphere’s magnetic dipole resonance, and at least one interfering eigenmode in 
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both the E' and E" responses is dominated by magnetic dipoles. In other words, we have 
eigenmodes in both irreducible representations that are dependent on a single parameter; 
the magnetic dipole polarizability of a single MgO-TiOo sphere. However, more generally, 
the correlation between a resonance of the single particle and modal interference is a recur¬ 
ring feature of the experiment. The single particle’s electric resonance at 6.5 GHz coincides 
with destructive interference in the E' response, between the |EU) and EW) eigenmodes. In 
Fig. ||)3 an d&, the single sphere’s magnetic quadrupole resonance at 7 GHz is also associated 
with a Fano resonance feature. It is interesting that the eigenmodes in our E" response ex¬ 
hibit, anomalistic, resonant behavior in the vicinity of the (neglected) magnetic quadrupole 
resonance. Indeed it is very likely that this frequency range should have significant mut- 
lipolar coupling channels [23j that have been omitted in our hybridization. However, for the 
purposes of the work here, the key result is the largest Fano resonance at 4.8 GHz, which 
is fully described by our hybridization theory. This Fano resonance is a realization of the 
predicted propensity that all-dielectric trimers have towards Fano resonances. 


IV. CONCLUSIONS 

We have presented an explicit study on the hybridization of optical responses in both 
plasmonic and all-dielectric trimers. We presented a simplified hybridization model to allow 
us to derive the eigenmodes of these structures analytically and observed the formation of 
modal interference in eigenmodes excited by normally-incident excitation. In this regard, 
plasmonic trimers were shown to exhibit nonorthogonality from retarded coupling channels, 
and dielectric trimers could further utilize bianisotropic coupling channels. A key prediction 
of our hybridization theory was then demonstrated experimentally: an all-dielectric trimer 
was shown to exhibit sharp, polarization-independent Fano resonances. The measurements 
were in good agreement with our simplified hybridization model, to therefore validate our 
approach. Our conclusions subsequently demonstrate that the use of full plasmonic hy¬ 
bridization is not necessary to deduce the dominant eigenmodes of multiparticle geometries. 
Yet, importantly, our analysis has demonstrated that hybridized eigenmodes can only be 
related to the set of intercoupled basis vectors in a simple form if the basis vectors are, 
in isolation, eigenmodes of their associated subsystem. Indeed, this raises a more general 
point of contention against many empirical applications of hybridization concepts and es- 
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tablishes why a more considered utilization of hybridization concepts in terms of subsystem 
eigenmodes is necessary. 
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Appendix A: Plasmonic trimer 


The dipole model we use for a plasmonic trimer assumes only electric dipoles from the 
individual nanoparticles. 


Pi = a B e 0 E o (ri) + a E k 2 ^ G 0 (r;, r ,-) ■ p ? 

iA 


(Al) 


Here, a E is the electric dipole polarizability of a particle and G 0 is the free space dyadic 
Green’s function, which acts on dipole moments as: 


G 0 (r',r) ■ p = 


0 ikR 


4t tR 


1 * 1 \ / 3i 3 

+ ’ ~ ’—IP ( + kR k 2 R 2 


kR k 2 R 2 ) 


(n ■ pjn 


where k is the wavenumber, n is the unit vector pointing from r to r' and R — |r — r'|. 
The associated eigenmode equation for a plasmonic trimer in this dipole model can then be 
written as: 


Vj = a E e 0 \vi + ^ a E k 2 G 0 (ri , rq) ■ v, (A2) 

iA 

If we define our eigenmode using the dimer and single particle eigenmodes as basis vectors, 
we can write any associated eigenmode as: 


I — a,i |p x ) + bi |B 3u ) + Cj |B lg ) + di |p y ) + ej |B 2u ) + /* |A g ) (A3) 

The six basis vectors are the dimer and single particle eigenmodes we define in the main 
text. The eigenmode equation, Eq. |A2[ can then be written as a 6 x 6 matrix equation in 
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terms of these basis vectors and the associated coupling channels. 
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where the coupling channels, A-D, are 

defined 

in Eq. [2] of the main text and Ab 3u , 


(A4) 


and Aa are the eigenvalues of the corresponding eigenmodes in the isolated dimer. 


AkR 
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(AS) 
(A6) 

(AT) 

(AS) 


These eigenvalues account for self-interaction of the dimer eigenmodes in the above matrix 
equation. The six solutions of Eq. |A4| provide the eigenmodes and eigenvalues for the 
plasmonic trimer in the main text. However, in regard to the nonorthogonality of |u 2 ) 
and 1 ^ 3 ), it is worth acknowledging that this nonorthogonality requires retarded coupling 


between the particles. Moreover, if refer to the matrix in Eq. A4 and the definitions of 
coupling channels in Eq. [2j we note that by neglecting retardation in the coupling between 
dipoles, we make the phase acquired through all coupling channels simply equal to that 


provided by the dipole polarizability. In other words, —, —, — and — all become real 
numbers. This then makes the matrix in Eq. |A4| proportional to a real symmetric matrix, 
and hence Hermitian, which then guarantees it has orthogonal eigenvectors. As we already 
have orthonormal basis vectors, this then makes all eigenmodes orthogonal. As such this 
shows that the retardation of coupling between particles is necessary for nonorthogonal 
eigenmodes. 
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Appendix B: Dielectric trimer 


For modeling the dielectric trimer, we use the coupled electric and magnetic dipole ap¬ 
proximation, which is described by the following two (coupled) equations: 


Pi = oi E e 0 E 0 (ri) + a E k 2 ( ^G 0 (ri, r,) • Pj - — V x Gofo, g) ■ nr, j 

V j^i C ° ) 

m,; = a//Ho(rj) + a H k 2 rj) ■ uij + c 0 V x Gofa, r,) ■ p ; 

Wi J 


(Bla) 

(Bib) 


where p, (m,) is the electric (magnetic) dipole moment of the z th particle, G 0 (ri,r,,) is 
the free space dyadic Greens function between the i th and j th dipole, a E (eta) is the elec¬ 
tric (magnetic) polarizability of a particle, Co is the speed of light and k is the free-space 
wavenumber. The extra bianisotropic coupling terms are given according to: 

AkR 


V x G 0 (r',r) • p = 


47ri? 


1 + M |nXp 


(B2) 


An eigenmode of the dielectric trimer, having electric dipoles v and magnetic dipoles u, will 


therefore satisfy the coupled electric and magnetic dipole model (Eq. Bl) as: 

Vj = a E e 0 Avj + a E k 2 ( ^G 0 (ri, rq) • v.,- - — V x G 0 ( r i; r,-) ■ J 

V C ° ) 

Ui = a H Xui + a H k 2 ^^G 0 (r*, rq) • uj + c 0 V x G 0 (rj, iq) • 


(B3a) 

(B3b) 


We can then use the dimer and single particle basis vectors from Fig. [3] in the main text to 
write an expression that will describe the eigenmodes of a dielectric trimer that are invariant 
under the dimer’s reflection symmetry operation. 


Wi ) — |p y ) + |B 2u ) + c' |A g ) + d[ |E' Z ) 


(B4) 


Analogous to the case for plasmonic trimers in 
4x4 matrix equation. 


Appendix [Aj we can rewrite Eq. B3 
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where the coupling channels, C-F, are defined in Eqs. [2] and [9] of the main text and the 
eigenvalue of the dimer eigenmode with ^-oriented magnetic dipole moments is: 


Ae' - ol h 1 + 


e ikR ( 2 ik 1 \ 

k + - — I 


47T.R 


K 2 J 


(B6) 


We can then reduce Eq. B5 to a 3 x 3 matrix equation when searching for only the doubly- 
degenerate eigenmodes given the eigenmode with radially-oriented dipole moments, |u 4 ), 


remains an eigenmode of Eq. B5 This is done by considering the eigenspace that is orthog- 


onal to u 4 ). In effect, this means we substitute b 
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= y/2 a' — \/3c' into Eq 

reduced matrix equation: 







( 1-V2C 

V3C-D 

-e 2 ^ 


f“ : ) 


f a B e 0 a' ^ 


-D 

oAa s 

0 


Cj 

= Aj 

a B e 0 c' 


y - Ei - \/2Fi 

V3Fi 

Q^Ae' j 




^ OiHd'i j 


(B7) 


From this expression, it is straightforward to find numerical solutions to this matrix equation 
and obtain both the eigenmodes and eigenvalues. However, as an addendum, we can consider 
the relations between the a', V and d coefficients of eigenmodes analytically. This will allow 
us to illustrate the interdependences of the collective eigenmodes on the resonances of the 


individual consitituent subsystems. Moreover, we can use the second row of Eq. B7 to first 
write a relationship between a' and d coefficients: 

-Da[ = (A i - A Ag )o!£eoc- (B8) 


From this relation, and using iq = \f2E L , the third row of Eq. B7 can also provide a 
relationship between d and d! coefficients: 

(Aj — Ae 'jDand'i = E\ (3(Aj — AA g )a,Eeo + Dy/6)d i (B9) 

These two relations are not sufficient to define eigenmodes analytically, because we need 


the first row of Eq. B7 to solve for Aj, which becomes quite nontrivial. However, for our 


purposes, Eq. B8 and B9 are sufficient to define the profile of an eigenmode for its given 
eigenvalue. Indeed, these two relations nicely illustrate the how the hybridized eigenmodes 
become dependent on each basis vector’s resonant properties and all the coupling channels. 
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